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ABSTRACT 

Cosmological hydrodynamical simulations of galaxy evolution are increasingly able to 
produce realistic galaxies, but the largest hurdle remaining is in constructing subgrid 
models that accurately describe the behavior of stellar feedback. As an alternate way 
to test and calibrate such models, we propose to focus on the circumgalactic medium. 
To do so, we generate a suite of adaptive-mesh refinement (AMR) simulations for 
a Milky- Way-massed galaxy run to z = 0, systematically varying the feedback im- 
plementation. We then post-process the simulation data to compute the absorbing 
column density for a wide range of common atomic absorbers throughout the galactic 
halo, including H I, Mg II, Si II, Si III, Si IV, C IV, N V, O VI, and O VII. The radial 
profiles of these atomic column densities are compared against several quasar absorp- 
tion line studies, to determine if one feedback prescription is favored. We find that 
although our models match some of the observations (specifically those ions with lower 
ionization strengths), it is particularly difficult to match O VI observations. There is 
some indication that the models with increased feedback intensity are better matches. 
We demonstrate that sufficient metals exist in these halos to reproduce the observed 
column density distribution in principle, but the simulated circumgalactic medium 
lacks significant multiphase substructure and is generally too hot. Furthermore, we 
demonstrate the failings of inflow-only models (without energetic feedback) at popu- 
lating the CGM with adequate metals to match observations even in the presence of 
multiphase structure. Additionally, we briefly investigate the evolution of the CGM 
from z = 3 to present. Overall, we find that quasar absorption line observations of the 
gas around galaxies provide a new and important constraint on feedback models. 

Key words: galaxies: formation - galaxies: evolution - methods: numerical - hydro- 
dynamics 



1 INTRODUCTION 

Following the flow of baryons in and out of galaxies is nec- 
essary to understand the star formation rate, stellar con- 
tent and overall structure of galaxies of all sizes. Significant 
progress has been made in our theoretical comprehension of 
the inflow side of this equation, with the improved under- 



standing of cold and hot-mode accretion (e.g. Binney 


1977 


Birnboim & Dekel|2003||Keres et al.|2005| Ocvirk et al. 


2008 


Kereset al. 2009 


Fauchcr-Gigucre et al.|2011 Van De 


Voort 


et al.||2011|). In this picture, halos with masses below about 



3 x 10 Mq accrete most of their gas in the form of ("cold") 



E-mail:chummels@astro. columbia.edu 



~ 10 K filamentary flows, while those that are more mas- 
sive develop an accretion shock and tend to accrete gas in a 
smooth, hot form. This implies that Milky Way progenitor 
systems, at high redshift, received baryons in the cold mode, 
but that present-day large spirals do not. 

Although the existence of these cold streams is robustly 
predicted by ACDM simulations, observational confirmation 
remains challenging. Using high-resolution SPH simulations, 



Faucher-Giguere fc Kere§| ( |2011[ ) demonstrated that the H 
I covering fraction of cold streams near damped Ly-alpha 
absorbers (DLAs) and Lyman-limit systems (LLSs) is 3% 
to 10% respectively at z ~ 2, an observational signal easily 
overwhelmed by galactic winds. Kimm et al.| ( |2011[ ) found 
that the predicted covering fraction of dense gas Nm > 10 20 
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cm" at z ~ 2 — 3 that would give rise to observable C II ab- 
sorption in background quasi-stellar objects (QSOs) is only 
about 5% for Milky- Way-mass halos. Despite these low cov- 



ering fractions, Pumagalli et al. (20111 found that the cold 



streams at high-redshift contribute significantly to the ob- 



served population of Lyman-limit systems (see also |Ribaudo| 
et al.|[20lT |, and that the transition from low to high mass 



should be accompanied by a disappearance of Lyman-limit 
absorbers at moderate to high redshift ( Ste wart et al.|201Tj ) . 
At low redshift, the inflow picture is even less clear, as accre- 
tion is predicted to transition out of the cold-mode. Detailed 
predictions are currently lacking, although see Joun g et al.| 
(20121; Fernandez et al. (20121 for a detailed comparison of 



H I gas properties between simulations and observations at 
z ~ 0. 

Galactic winds are the other side of the baryon flow 
equation ( Dave et al.|2012 l, and evidence from the current 
stellar content of galaxies indicates that these must carry a 



substantial amount of mass and metals ( Mandelbaum et al. 
20(36} |Oppenheimer fc Dave|2008||Behroozi et al.|2010[ ). Nu- 
merical simulations have struggled to reproduce these large 



mass outflows (Katz 1992 Abadi et al. 2003 Hummels & 



Bryan 2012J, but the addition of subgrid models can pro- 



duce high mass winds (e.g. Springel fc Hernquist |2003[ | Stin- 
son et al. | |2006l lOppenheimer & Dave 2006; Scannapieco 



et al.|2006| |Governato et al.|2012 l. However, a detailed un- 



derstanding of what drives winds is still lacking, with even 
the basic physical mechanism (e.g. ram pressure from super- 
nova bubbles or radiation pressure) in dispute (e.g. |Murray| 



et al.|2011||Sharma fe Nath|2012||Hopkins et al.|2012[ ). 

One important way to probe both the inflow and out- 
flow of baryon, and an important test for any cosmological 
galaxy formation model, is the observed correlation between 
absorption lines in background quasars and the projected 
distance (and velocity offset) to the associated galaxy. A 
substantial body of such data has been built up at high red- 



shift (Simcoe et al. 2004 Steidel et al. 20101. These data 



are beginning to be compared against simulations at high 



redshift (e.g. |Altay et al. 


2011 Tescari et al. 


2011 


Fuma- 


galli et al.|2011||Barnes et al.|2011||Goerdt et al.|201 


2||Shen 


et al. 2012|). At low redshift, where we focus, there exists 



a large body of data combining galaxy imaging and QSO 



spectroscopy dChen et al. 2001 
& Chen||2008| |Yao et al.||2008 ' 


|Chen fe Tinker|2008[|Thom| 
|Chen et al.||2010| |Prochaska 


et al. 


2011 


Thorn et al.||2011 


Tumlinson et al.||2011 1, but 



only a few predictions from numerical simulations ( |Stinson| 
et al.|2011||Fernandez et al.|2012 1. 



In this paper, we use high-resolution cosmological sim- 
ulations of a galaxy run to z — to make predictions for a 
wide variety of absorption lines as a function of impact pa- 
rameter and compare them against observational samples. 
We include tracers of cold gas (T ~ 10 4 — 5 x 10 4 K) such as 
HI, Mg II, Si II, Si III, and Si IV, probes of warm-hot gas (C 
IV, N V, and O VI, for T ~ 10 5 - 3 x 10 s K), as well as gas 
at the virial temperature of the halo (O VII; T ~ 10 6 K). 
We will show that these comparisons make for a sensitive 
(and challenging!) probe of the phase of the inflowing and 
outflowing gas. These runs are conducted on identical initial 
conditions, systematically varying the stellar feedback pre- 
scription to demonstrate their effect on the column density 
distribution of absorbers in the galactic halo. In Section [2j 
we describe the details of our simulations and our post- 



processing steps on both the numerical and observational 
datasets. Section [3] presents the results of these simulations 
and the consequent analyses. We discuss these results and 
their implications for the nature of outflows, feedback and 
for use as a probe of different feedback models in Section 
[4] Lastly, we summarize these conclusions in Section [5] and 
discuss future directions for subgrid models to match these 
added observational constraints. 



2 METHODOLOGY 
2.1 Simulation Code 

The numerical methodology and initial conditions used in 
this paper are very similar to those employed in |Hummels| 
& Bryan (2012 hereafter Paper I), and so we only briefly 



describe them here. We perform all of our simulations us- 
ing the adap tive mesh refinement (AMR) hydrodynamic s 
code, EnzcQ ( |Bryan fc Norman|[l997| |Q'Shea et aL]|2004 |. 
Enzo uses AMR to dynamically restructure the grid over 
the course of the simulation to achieve higher resolution 
in regions where predetermined refinement criteria are met. 
It treats dark matter and stellar populations as collision- 
less particles modeled with an N-body particle-mesh gravity 
solver ( Efstathiou et al.[T9 85 ) . The equations of gas dynam- 



ics are solved using the ZEUS hydrodynamics engine ( Stone 
fc Norman|T992"l ) 



A minimum pressure floor is included in these simula- 
tions in order to keep the Jeans length resolved across all 

We follow the implemen- 
(|2001|), wherein the ratio of the 
Lj /Ax > 8 by 



grid cells (Truelove et al. 



tation of Machacck et al. 



Jeans length to cell size is maintained at J : 
the addition of an artificial pressure. 

We model the process of star formation following a pre- 



scription similar to that in Cen & Ostriker ( 1992 1. A cell will 



spawn a star particle if the following criteria are met: (1) the 
cell's gas overdensity exceeds 5sf, (2) the local velocity di- 
vergence is negative, and (3) the cooling time is less than 
the dynamical time for the gas cell. When these conditions 
exist, a star particle with mass m* is created according to 
the relation: 



At . 3 
m, = e SF - PgasAx 

Edyn 



(i) 



where esf is the star formation efficiency per dynamical 
time; At is the local timestep; tdyn is the local dynamical 
time, defined as td yn = (37r/32Gp) 1,/2 ; p gas is the gas den- 
sity; and Aa; is the grid cell width. To prevent a large num- 
ber of low mass particles, we adopt a minimum star particle 
mass of m*. m i n . If m» < ra*, m i n , the star is formed with a 
probability of m*/m*. m i n , and its resulting mass is the min- 
imum of m^min and 90% of the mass in the gas grid element 
(e.g. Tasker & Bryan 20061. For these simulations, we set 



m*, mi n = 10 s M . We find that e S F = 10~ 2 and <5 SF = lO 3 
approximately reproduces the Kennicutt-Schmidt Law (see 
Paper I). This value of Ssf is equivalent to 0.1 cm -3 at z = 3. 



http:/ /enzo-project.org 
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2.1.1 Supernova Feedback 

Our implementation of stellar feedback includes a prescrip- 
tion for prompt injection of thermal energy from type II 
supernovae to the surrounding medium. Individual star for- 
mation events produce star particles of m* ~ 10 4-5 Mq, 
each representing an entire stellar population. The high- 
mass stars within this stellar population will become su- 
pernovae and send energy, mass and metals into the local 
medium. Because a stellar population is rarely created in- 
stantaneously in nature, we "activate" a small component 
of the star particle's mass each timestep to represent the 
stars that were effectively created during that period. We 
parameterize the extended timescale over which this occurs 
with: 



(2) 



At t — to _( t _ to )/ T 
ttifb = m, — e 

T T 

where ttifb is the mass of newly activated stars for a 
given timestep, to is the star particle creation time, and 
r = max(tdyn, 10 Myr). In order to prevent unrealistically 
short dynamical times, we apply a floor of 10 Myr to r. We 
normalize our activated star mass so that the total mass re- 
leased over the lifetime of the star particle is one-fourth the 
star's original mass (i.e. ^7TIfb = 0.25m*). 

The returned energy eFB and returned metals m me taia 
due to stellar feedback are directly proportional to the mass 
of these new stars: 



eFB = 771FBC 6FB 



etals = "iFB V 



(3) 
(4) 



where c is the speed of light; efb is a feedback efficiency 
parameter; and y is the fraction of metals returned relative 
to the activated stellar mass. Thus at any timestep At af- 
ter the creation of a star particle, it returns jtifb mass, eFB 
energy, and m mcta i s metals to the cell containing the star 
particle. This thermal energy should heat the surrounding 
gas, increasing the Jeans mass, and prevent further imme- 
diate star formation. 

In the present investigation, we use a value of y = 0.02 



consistent with previous work (Smith et al. 20111, and we 



include metal and mass feedback in all of our simulations. In 
four of our simulations, we vary the thermal energy efficiency 
value: e F B = 3x 10~ 6 (low), lx 10" 5 (medium), and 3x 10 -5 
(high). These different energy efficiencies correspond to one 
10 51 erg SN for every 180, 60, and 18 M stars formed re- 
spectively, effectively sampling the range of feedback effi- 
ciencies used by other studies on this topic (see discussions 
in Section 4.3 and Paper I). 



2.1.2 Cooling Suppression 

One of the shortcomings of the thermal feedback model is 
that thermal energy dumped into an underresolved volume 
will fail to bring that gas up to a high enough temperatures 
to prevent it from immediately radiating away this energy. 
To address this, many researchers have adopted an ad hoc 
solution, to turn off radiative cooling in the gas for a pe- 
riod after a star formation event ( Gerritsenj|1997| |Thacker 

& Couchman 2000 ; Sommer-Larsen et al. 2003 ; Stinson et al. 

^ 

2006; Gover nato et al.|2007[ ) This model is justified as an ap- 
plication of the Sedov- Taylor solution for a supernova blast 



wave applied at scales below the spatial resolution of the 
simulation, although these arguments are not entirely phys- 
ically consistent (see Paper I). However, cooling suppression 
seems to effectively quench star formation and slow the over- 
cooling problem particularly when combined with an addi- 
tional feedback prescription ( Colm et al.|2010| Agertz et al.| 
[20TT] |Guedes et aLl|20"TT| |Piontek fc Steinmetz||2011| Pa- 
per I). 

Our implementation of cooling suppression is relatively 
simple. After a star-formation event, we shut off cooling for 
a period t SU p P in the gas cell containing our newly-formed 
star particle. We allow this effect to track the star particle 
until taupp has passed, just as our thermal prescription for 
feedback tracks the star particle for a time after its creation. 
In this study, we set t a u PP = 50 Myr, a value consistent 
with previous studies ( Stinson et al.|2006 Colm et al.|2010 
Paper I). 



2.1.3 Radiative Heating and Cooling 

In modeling the thermal state of the gas in these simulations, 
Enzo accounts for photoelectric heating from a metagalac- 
tic ultraviolet background and radiative cooling from metals 
present in the gas. Because directly modeling the fractional 
components in each ionic state of every element is too com- 
putationally costly for simulations, we adopt a different tac- 



tic. We employ the method pioneered by Smith et al. ( 2008 1 



and subsequently used by other groups in a variety of con- 
texts (Booth & Schaye 2009), wherein all of the phases of H 
and He are followed individually throughout the simulation, 
but for more massive elements, a multi-dimensional lookup 
table is consulted to determine their aggregate contribution 
to heating and cooling using a single all-encompassing metal- 
licity field. All phases of H and He (H I, H II, He I, He II, 
He III, and e~) are calculated by a non-equilibrium, primor- 
dial chemistry network, which uses a backwards-differencing 
formula to account for collisional ionization, recombination, 
brehmstrallung, Compton cooling, photoexcitation and pho- 
toionization rates in these species ( |Abel et al.|1997| |Anninos| 
et al.|19 97). The metal heating/cooling lookup table is com- 
piled usin g version 07.02.01 of the spectral synthesis code, 
Cloud^ jFerland et al. 19981. The table consists of five 
dimensions, one each for density, temperature, metallicity, 
electron fraction, and redshift. The redshift index is included 
to determine the intensity and spectrum of the metagalactic 
ultraviolet background at that time, which we activate at 
z = 7 and henceforth follow the prescription set in IHaardt] 
fe Madau] ( |1996[ ) as updated in this version of Cloudy. Prior 
to reionization at z = 7, metal heating and cooling are only 
due to collisional ionization. 



2.2 Initial Conditions 



In this study, we adopt the WM AP 5-year results ( Komatsu 
et al.||2009| , specifically adopting Q = 0.258, ft a = 0.742, 
ftbaryon = 0.044, cr 8 = 0.796, H a = 71.9 km s" 1 . Using the 
mits program included in the Enzo distribution, we generate 
a 128 3 particle grid over our 20h~ 1 Mpc volume and carry 
out a low-resolution simulation to identify a dark matter 

2 http://www.nublado.org 
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Name 


Description 


M 20 o 


Mum 


Mstar 


M gas [hot, cold] 1 


^200 


^circ,max 






^supp 






no 10 MrO 


no 10 m^i 


MO 10 M^"l 


MO 10 M^l 


fkncl 
l K Pc; 


( km s ) 








V 


Vanilla 


110 


93 


17 


2.5 [1.8, 0.6] 


210 


590 










LP 


Low Feedback 


110 


92 


15 


3.6 [2.8, 0.8] 


210 


560 


3 x 10~ 


-6 





MF 


Medium Feedback 


110 


91 


11 


5.1 [4.4, 0.8] 


210 


520 


1 x 10~ 


-5 





HF 


High Feedback 


90 


84 


2 


4.3 [3.0, 1.3] 


200 


200 


3 x 10- 


-5 





MFCS 


Medium Feedback & 


100 


90 


7 


6.6 [5.7, 0.9] 


200 


390 


1 x 10- 


-5 


50 


Cooling Suppression 







Table 1. We present each of our simulations, their simulation parameters, and the resulting halo characteristics at z = 0. 



halo with mass M200 ~ 10 12 Mq (see Paper I for more details 
- we use halo 26 of that paper). Once the halo is identified, 
we select all particles within 2 virial radii at z = 0, and trace 
them back to their original positions at z = 99. The region 
containing all of these particles at z = 99 is further refined 
with two nested refinement regions, each with 2-times better 
spatial resolution and 8-times lower mass dark matter parti- 
cles, resulting in Mum = 4.9 x 10 6 Mq dark matter particles 
in the vicinity of our target halo. Additional (adaptive) re- 
finement occurs in a cell when its gas mass (or dark matter 
mass) exceeds 4 x 10 6 M Q (or 1.9 x 10 7 M for the dark mat- 
ter), resulting in a further increase of the spatial resolution 
by a factor of two, up to nine levels of refinement, achieving 
a maximum resolution of 425 co-moving pc. The timescales 
over which our particles of different mass (i.e. stellar and 
dark matter) suffer from relaxation effects is on the order of 
a Hubble time, assuring these spurious effects are avoided. 



2.3 Description of the Simulations 

We performed five different simulations using identical ini- 
tial conditions and simulation parameters as described 
above. The way in which these simulations differed was sim- 
ply in their supernova feedback and cooling suppression pa- 
rameters. These runs include: a control run with no energetic 
feedback or cooling suppression ["Vanilla"]; three runs with 
no cooling suppression but varying feedback efficiencies of 
e FB = 3 x 10" 6 , 1 x 10~ 6 , and 3 x 10" 5 ["Low Feedback", 
"Medium Feedback", and "High Feedback", respectively]; 
and a run including cooling suppression and feedback with 
£fb = 1 x 10 -6 ["Medium Feedback and Cooling Suppres- 
sion"]. These runs and some of the properties of their halos 
at z — are listed in Table [l] The simulations produced 
in this study have identical initial conditions and physical 
parameters to those associated with run H26SPM (i.e. the 
canonical run) in Paper I save for the treatment of feedback 
and the presence of metals. See that work for more detail. 
It is worth mentioning that we do not tune our models in 
any way to produce the results herein, all parameters be- 
ing either deduced from theoretical calculations or values in 
the literature. Unless otherwise noted, all units herein are 
in comoving units, where h -1 is excluded. 



2.4 Analysis Code 

We applied a series of post-processing analyses to each sim- 
ulation so as to render synthetic observations comparable 



to observational datasets for these simulated galaxies. Each 
simulation outputs a datafile to disk every 10 Myr, repre- 
senting the state of the gas, stars and dark matter at that 
time. Following the method described in ( |Smith et al.|201l| ), 
we processed each datafile to compute the spatial distribu- 
tion of various different ionic species. The density of the ith 
ionic species of a given element, X, is given by 



nx z = n H x (nx/nu) x (n x Jn x ) 



(5) 



where hh is the density of hydrogen, (nx/nn) is the ele- 
mental abundance of X and (nxi/nx) is the ion fraction 
in the ith energized state. We internally track the neutral 
and excited states of hydrogen and helium throughout the 
simulation, but we represent all metals as a single metal- 
licity field. Thus, we assume the metallicity of any given 
metal X tracks with this overall metallicity Z consistent 
with solar abundance ratios: nx/nu = Z x (nx/nii)©. To 
calculate (nx t /nx), we create a lookup table using the spec- 
tral synthesis code, Cloudy (version 07.02.01), which gives 
(nxi/nx) uniquely determined by the gas temperature, gas 
density, and the photoionizing spectrum and intensity at 
that time. We use the redshift-dependent ultraviolet back- 
ground included in Cloudy as the photoionizing source (a 



2005 updated version of Haardt & Madau 1996 1, which takes 



into account reionization by QSOs and galaxies. The Cloudy 
lookup table accounts for collisional ionization as well as 
photoionization by the metagalactic UV field (although we 
do not track local photoionizing sources). This infrastruc- 
ture coupled with the yt simulation analysis suit^] ([Turk 
et al.|20lT |, allows us to add any arbitrary ion and element 



to our simulation as an additional data field, as defined by 
the local baryon density, metallicity, temperature, and red- 
shift. We process all of our data outputs to include fields for 
the astrophysically important ions of H I, Mg II, Si II, Si 
III, Si IV, C IV, N V, O VI, and O VII. 

For each data output of each simulation, we locate the 
galaxy of interest, and create a series of ten column-density 
projections centered on the galaxy from randomly oriented 
vectors passing through the galaxy in each of the selected 
fields. Each projection consists of two 128 2 -pixel images, one 
with a field of view of 2 Mpc and the other with 100 kpc, al- 
lowing us to achieve a large dynamical range in resolution (1 
Mpc to 1 kpc resolution in the central region of the galaxy). 
These images provide a map of the column density of each 
species and can be compared directly against observations. 

3 http://yt-project.org 
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Figure [T] displays a few such projections from the MFCS run 
in H I, C IV, VI and total gas density at z — 0.5. For il- 
lustration's sake, these images are 500 kpc on a side, with a 
circle overplotted representing r2oo, the radius at which the 
galactic density drops to 200 times the mean density of the 
universe. The column density limits chosen are those used 
later in the radial profiles for each species. 

For the image plane, pixels are identified by their co- 
moving distance from the central pixel. These pixel radii are 
then binned, and the bins are populated by the column den- 
sity maps of a series of images for a given galaxy. The median 
column density value in each radius bin is computed to pro- 
duce a radial profile for that galaxy, representing the typical 
column density one might expect to detect at a given radius 
were one to observe it from a random angle. We found early 
on that there were significant temporal peculiarities (e.g. in- 
falling satellites) which made it difficult to compare one run 
at a given redshift against another, so we made radial pro- 
files consisting of data outputs over a range of redshifts to 
wash out this effect. To give some insight into the range of 
column densities at each radius, we also include some quar- 
tile column density ranges (i.e. 25% and 75%) in these radial 
profiles. 

For the physical quantities of metallicity and temper- 
ature, we take a mass-weighted average projection for each 
quantity as opposed to a column density. All further process- 
ing steps on these projections follow the same prescription 
as for the atomic species projections. 



2.5 Processing the Observational Datasets 

In order to say something meaningful about our numeri- 
cal models, we compare them directly against observations 
taken in a similar manner. Typically, gas in these species is 
not observed in emission, but rather in absorption due to its 
low density (emission oc n 2 , whereas absorption oc n). We 
identified observational studies which correlate the proxim- 
ity (i.e. p = impact parameter) of absorbers to nearby host 
galaxies for a subsample of our chosen atomic species. The 
observational studies against which we compare our data 



are: H I (Lya) 


from |Prochaska et al. (20111 (utilizing data 


from|Penton et al. 2004; Danforth & Shull 2008; Tripp et al. 


2008); Mg II: 


Chen et al. 


( 


2010 


1; C IV from Chen et al. 



Yao et al. 



J2008| ). 

In two of these studies (e.g. Chen et al.|2001 20101, the 
absorption features were not fully resolved into individual 
components, and so only equivalent widths of the absorbing 
gas were provided, whereas in all of the other studies both 
equivalent width and column densities of the absorber are 
reported. In order to directly compare our results against 
all of these studies, we needed to convert any equivalent 
widths to column densities. Converting between equivalent 
width and column density is not a trivial proposition, as it 
requires knowledge of the temperature and velocity distri- 
bution of the gas containing the absorbing species. To make 
an approximate translation, we used a "curve of growth" 
analysis. We applied a fitting function relating equivalent 



W 



c 1+t /(2v/2) 
(t) 2 -(lS 

2 2 Nifi u \i u 



+ 



b 7A (tq -1.254) 



1/2 



tq < 1.254 
to > 1.254 



(6) 



(7) 



where W is the equivalent width of a given absorption line 
for a transition from state I to state u, b = is the 

velocity dispersion of the absorbing gas, c is the speed of 
light; ro is the optical depth at absorption line center, 7 is 
the intrinsic width of the absorption line, A is the wavelength 
of the transition, Ni is the column density of material in the 
I state, and fi u is the oscillator strength of the transition. 

This relation converts equivalent width to column den- 
sity, but we want the inverse, so we numerically inverted 



the function applying the bounds-bisection method (Press 



et al. 19921. We used Verner et al. (19961 for accurate os- 



cillator strength, wavelength, and 7 values for the various 
transitions. 

To deduce the appropriate value of b for our absorber 
populations, we utilized the equivalent width and column 



density data in Tumlinson et al. ( 2011 1 and Prochaska et al. 
( |2011 1. Applying equation |6j to calculate column density 
from these equivalent widths for a variety of fr-values, we 
then made a x-squared fit with the listed column density 
values to ascertain the best choice of b for O VI and H I 
respectively. Our results indicated 60 vi = 58 km s^ 1 and 
&hi = 29 km s . Since the gas responsible for absorbing C 
IV is in a similar temperature regime (slightly lower) to that 
of O VI absorbing gas, we used feciv ~ i>ovi = 58 km s _1 . 
Additionally, Mg II has a similar ionization energy to H I, 
and Mg II is though t to trace H I at Nun > 10 16 cm -2 
iPutman et al. 20121, so we let &M g n ~ &m = 29 km s . 



width and column density as found in Draine (2011) 



Because the lines of these species may have multiple com- 
ponents, these b values should not be take too literally, but 
for the purposes of this current comparison, they suffice. 

The uncertainties endemic to this type of conversion 
are mostly contained in our choice of b value. Changes in 
b of 3% result in changes in calculated column density of 
~ 3% and ~ 25% depending on whether the absorber is in 
the linear or saturated region of the curve of growth respec- 
tively. Fortunately, the majority of our C IV (95%) and Mg 
II (70%) absorbers lie in the linear regime, with only the 
strongest absorbers populating the saturated region. It is 
worth noting that we are only applying this conversion to 
two observational datasets (Mg II and C IV) for compari- 
son's sake, and not to our own numerical datasets, so any 
problems with this conversion have limited repercussions to 
the study at hand. Furthermore, the main results of this pa- 
per come from comparisons against the O VI dataset, where 
this conversion was not required. We thus convert the ob- 
servational datasets to approximate column densities (and 
proper units to comoving units) for comparison against the 
numerical results. 



3 RESULTS 

Our suite of simulations yields halos with similar bulk prop- 
erties, consistent with the fact that they all share the same 
initial conditions. We catalog the mass of each component 
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Species 



Redshift 



Reference 



H I 

Mg II 
Si II 
Si III 
Si IV 
C IV 
N V 
O VI 
O VII 



13.6 eV 

15.0 

16.3 

33.5 

45.1 

64.5 

97.9 

138. 

739. 



0.05 - 0.2 
0.1 - 0.5 



0.08 - 0.9 

0.1 - 0.35 




Prochaska et al. 



I 2011 1 



Chen et al. " pOlO I 



Chen et al. 



(2001) 



Tumlinson et al. 



120111 

v;u. et ai. ^ns^r 



Table 2. The observational datasets used to compare against our 
numerical radial profiles. We also include ionization energies for 
our selected ionic species. 



(dark matter, gas, and stars) and the maximum circular ve- 
locity for the primary halo in each of our simulations at 
z = in Table [Tj The mass in each component is measured 
interior to r2oo at z = 0. There is little variation in the dark 
matter and total mass across our simulations except in the 
case of the high feedback run where large baryonic outflows 
have softened the potential and diminished the overall mass 
of the system. 

In Paper I, we examined the rotation curves and star 
formation histories of a set of simulations very similar to 
those used in this paper. Here, we focus on the observable 
properties of the circumgalactic medium (CGM) at low red- 
shift - in particular, we report on the column density dis- 
tributions of a wide variety of atomic species, observable in 
absorption against background quasars. We initially restrict 
ourselves to low redshift, as there exists a substantial body 
of observations which link absorption strength to projected 
distance from the galaxy. In the following sections, we first 
examine the column density distributions, then explore the 
physical properties of the gas giving rise to this absorption, 
and finally investigate the redshift evolution of the predicted 
absorption strength. 



3.1 Radial Profiles for Atomic Species 

In this study, we examine the radial column density distri- 
bution for a number of absorption-line-generating species, 
well-sampling the range of ionization energies (and thus tem- 
perature/density regimes) of the CGM as shown in Table [5] 
At the two ends of the spectrum, the H I traces the cold, 
dense clouds, whereas the O VII probes the hot coronal gas. 
Between these regimes is a continuum of layers of mate- 
rial acting as the interface between these two phases. This 
layering effect takes place on many different spatial scales 
from small star-forming clumps of cold gas up to the galactic 
scale. Figure [I] demonstrates this at the latter scale, where 
the cold H I of the inner halo is enshrouded in warmer C 
IV-bearing gas, which is further encased and extended by 
gas exhibiting O VI. This provides us with some intuition 
for understanding the following galactic profiles for these 
species. 

We produce radial profiles of the median gas column 
density in nine different ionic species as a function of pro- 
jected radius extending from 2 kpc to 1 Mpc for our simu- 
lated galaxies. We probe this gas from the disk-dominated 



region (r < ri nncr ~ 20 kpc) through the CGM (ri, mcr < 
r < r2oo) into the intergalactic medium (r > r2oo)- Figures 
[5] [3] and [4] present these results (Figure [5] shows total gas 
column densities for comparison). The left and right sides 
of each figure show the radial profiles of the same quantity 
but for different simulations. The left sides display a com- 
parison of simulations with and without different forms of 
feedback: V, MF and MFCS, whereas the right sides exhibit 
a progression of thermal feedback intensities: LF, MF and 
HF (the MF simulation is repeated on both sides to ease 
comparison). Each profile gives the median value expected 
in an observation of a galaxy at a given radius over the red- 
shift range z — 0.1 — 0.5 (selected to approximately match 
the redshift ranges of the observations). To give some idea of 
the variation around the median absorber strength, in each 
panel we overplot a color band representing the quartiles of 
the column density distribution (i.e. 25% to 75%) for one 
model: MFCS (left) and HF (right). The quartile distribu- 
tion associated with the HF model is most indicative of the 
spread of all the other models, whereas MFCS produces a 
particularly spread-out distribution suggesting more multi- 
phase substructure. In Figure [2] we present radial profiles for 
our lowest ionization energy species, H I, Mg II, and Si II; 
Figure [3] shows medium energy ions, Si III, Si IV, and C IV; 
and Figure [4] displays the highest energy species N V, O VI, 
and O VII. Overplotted on several profiles are low-redshift 
observational datasets, where magenta points represent de- 
tections and grey arrows mark limits. These datasets, their 
sources, and their redshift ranges are listed in Table [2] 

The H I, Mg II, and Si II datasets (the H I group) 
in Figure [2] demonstrate similar behavior across all of our 
models, which is consistent with the fact that their ioniza- 
tion energies are very close to each other (see Table j2|. All 
of the simulations show a high-column-density core (i.e. the 
disk) in these species out to r ~ 10 kpc, at which point the 
profiles drop to trace levels by 10-50 kpc (i.e. the extent of 
cool gas surrounding the disk). This transition radius, ri nncr , 
is related to the extent of the galactic wind in that it marks 
the outer boundary of this drop, where the material falls 
off to the ~ 10 13 cm -2 level. Each feedback model predicts 
a different n llnC r, sometimes differing over 40 kpc in radial 
location, identifying ri nncr as a potentially useful differen- 
tiator between feedback models. Our V and HF simulations 
can match the bulk of the Mg II observations; however none 
of the models reproduce the most extended high column 
density observations in H I and Mg II. Note that we have 
subsampled the H I observations ( Prochaska et al.|2011 1 in- 
cluding only those detections associated with galaxies in the 
0.5L* < L < 2L* range since these best match the charac- 
teristics of our simulated galaxy. 

In Figure [3] the medium-ionization species Si III, Si 
IV, and C IV (the "C IV group") display somewhat sim- 
ilar behavior to those of the H I group, each exhibiting a 
core and precipitous drop in material around 10-40 kpc. The 
floor level of the C IV absorption strength (beyond rinner) 
is larger than either of the other two lines, presumably due 
to the fact that it is tracing warmer gas. We also compare 
to the observed C IV dataset, but note that this compari- 
son is difficult due to the somewhat contradictory limits and 
detections in the outer halo; however, all of the simulations 
appear to be an order magnitude too low in column density 
in the probed radial regime (r ~ 10 — 200 kpc). 
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Figure 1. Sample projections from the medium feedback and cooling suppression run at z = 0.5. Each image represent a column-density 
map of a region with 500 comoving kpc on a side. The white circle denotes r-200 = 160 kpc. 



The tracers of warm and hot gas, N V, O VI, and O 
VII, have more extended and flatter distributions than the 
H I and C IV groups. The O VII distribution, which probes 
gas with temperature close to the virial temperature of the 
halo is particularly flat (note the different vertical scale). 
Our models that best fit these observational datasets are 
the HF and MFCS simulations; however even those profiles 
are an order of magnitude below the observations of O VI 
in the r ~ 10 — 100 kpc range. These same models reason- 
ably fit the O VII observations, although this dataset has 
relatively weak constraints. The vanilla run conspicuously 
falls off in column density at r ~ 50 kpc in all of the species 
(we will address the reason for this in more detail below). 
The O VI observations are separated into "passive" galaxies 
(plotted with an asterisk) and star-forming galaxies (plot- 
ted without an asterisk) (see |Tumlinson et ak]|2011[ for a 
full description) . While one would expect the active galaxies 
from the observational sample to be more like our simulated 
star-forming disk galaxy, it appears that the passive systems 
yield a much better match. 



3.2 Radial Profiles for Physical Quantities 

In order to better understand the physical basis for the sim- 
ulation predictions, we next turn to the physical properties 
of the CGM gas itself. We present radial profiles for median 



gas column density, median metallicity, and median temper- 
ature in Figure [5] in a manner similar to Figures [2] |3J and 
[4] In order to make a fairer comparison to the projected 
column density profiles, we compute these profiles using the 
same machinery, first carrying out several projections for 
each quantitiy (mass-weighting for the metallicity and tem- 
perature), and then computing the median value for each 
projected radius over the same z = 0.1 to 0.5 redshift range 
used earlier. The profiles computed in this way are similar 
to the three-dimensional quantities, but are more straight- 
forward to compare at a given projected radius. 

For the gas column density profile, there is not a great 
deal of variation (within a factor of 2) from one simulation 
to another, suggesting that different feedback models play 
only a modest role in rearranging the radial distributions of 
the gas mass within a galaxy. 

The metallicity profiles demonstrate a much more dra- 
matic difference between those runs incorporating energetic 
feedback (LF, MF, MFCS, HF) and the run which does not 
(V). While all of the models share roughly the same median 
metallicity in their inner 50 kpc (within 1 dex of Solar), 
the median metallicity of the vanilla run plunges to nearly 
primordial values beyond that distance, explaining why the 
vanilla run drops off in C IV, N V, O VI and O VII. All of 
our simulations share the ability to produce and return met- 
als to the surrounding medium from star-formation events, 
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Figure 2. Radial profiles for column densities of H I, Mg II, and Si II in comoving units covering the period z = 0.1 — 0.5. The color bands 
represent the quartile values (±25%) overplotted on the median for two of the models: MFCS and HF. Also shown are observational 
datasets with detections as magenta points and limits in gray arrows from the datasets listed in Table[2] The numerical datasets contained 
in Figures [2] [3] [4] and [5] are available in tabular form for public download at http://chummels.org/CGM.html . 
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Figure 3. Similar to Figure [2] but for Si III, Si IV, and C IV. 



but only the "feedback" runs are capable of depositing en- 
ergy along with those metals. It is clear that the inclusion of 
some form of energetic feedback is necessary to mix and dis- 
perse high-metallicity gas beyond 50 kpc. Lacking energetic 
feedback, the V simulation effectively has no outflow of ma- 
terial, and the metals get mixed up solely by gas motions in 



the inner halo region. This fact makes it very useful for des- 
ignating as an inflow-only model, thus cleanly isolating the 
characteristics of infalling material versus models generating 
outflows. Because our inflow-only simulation shows similar 
concentrations of gas in the outer halo, we can conclude that 
most of the gas in the halo must have been deposited there 
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Figure 4. Similar to Figure [5] for N V, O VI, and O VII. The O VI dataset (Tumlinson et al,|2011[ | is split into star-forming gala 
(circle symbols and plain arrows) and passive galaxies ("asterisk" symbols and "asterisk" arrows). 



by inflows. Evidently, the main effect of feedback is to mix in 
higher-metallicity material to large halo distances, but not 
to enhance the density there dramatically. 

All of the simulations share the same general behavior 
in the projected temperature profiles: low in the core (the 



cold gas of the disk) out to ~ ri nn er> higher values extend- 
ing to around the virial radius, and then decreasing beyond 
that. Note that the "cold" gas temperatures (T ~ 10 5 K) 
seen in the core in these profiles reflect the presence of multi- 
phase gas. The LF, MF and MFCS runs exhibit a high peak 
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Figure 5. Radial profiles for physical quantities: column density of gas, mass-weighted metallicity, and mass-weighted temperature 
covering the period z = 0.1 — 0.5. The color bands represent the quartile values (±25%) overplot on the median for two of the models: 
MFCS and HF. 



temperature of T ~ 20 million Kelvin at n nncr , where the 
high density of the galactic disk drops off. Conversely, the 
V and HG runs peak farther out (r ~ 40 kpc) at a lower 
temperature (T ~ 7 x 10 s K). The lack of clear trends in 
the temperature profiles indicate that low values in the cen- 



ter reflect simply how much cold gas survives in each model, 
while the temperature in the outer parts of the halo probably 
depends largely on the halo's virial temperature. These tem- 
perature profiles, however, are clearly related to the species 
profiles shown earlier - the low ionization species (particu- 
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larly the H I group) only appear where the mass-weighted 
gas temperature is low. Even the C IV group are correlated 
strongly with the extent of cold gas, probably indicating that 
it is formed in multi-phase gas. Finally, the column density 
profiles of the higher ionization species are generally flatter 
than the gas profile, indicating that they are primarily com- 
ing from the outer halo and their slight overdensity in the 
core is mostly a projection effect. This is particularly true 
of the O VII profile, which is very flat. 



3.3 Redshift Evolution of Radial Profiles 

There exists a significant body of observations as well as a 
growing number of computational studies investigating the 
radial profiles of metal-bearing gas in the CGM of galaxies 
at high redshift (e.g. |Altay et al.||2011| [Tescari et al.||2011 



Fumagalli et al.|2011| |Barnes et al.|20li| |Shen et al.|2012 | 
While the focus of the current work is low redshift, it is 
worth examining how the CGM evolves over time in our 
simulations. Figure [6] presents radial profiles of gas column 
density, mass-weighted metallicity, and O VI column density 
for two of our models (V and MFCS) as a function of red- 
shift. These profiles are constructed using the same method 
of median-sampling many data outputs over a range of red- 
shifts so as to wipe out any temporal peculiarities in the 
profile. The figure demonstrates trends in the behavior of 
certain quantities with time sampling redshifts from z = 4 
to z = 0. 

The gas column density behaves similarly in both of our 
models, consistent with an expanding universe. From an av- 
erage redshift of z = 3 to z = 0, the scale factor changes 
by a factor of four, so that the column density of a uniform 
medium integrated over a constant comoving path length 
should decrease by a factor V 2 = 16, which is approxi- 
mately consistent with these profiles, particularly at large 
projected radii. On the other hand, metallicity appears to 
increase over time in both models, especially for the inflow- 
only simulation, although the effect is quite weak. In the 
MFCS run, the metallicity grows with time at all radii, but 
by only a few percent. Much of the star formation occurs 
at early times and is spatially distributed throughout the 
galaxy, quickly enriching the gas at all radii, so that subse- 
quent star formation seems to do little to further increase 
the metallicity throughout the halo. 

Finally, the O VI radial profile in this figure exhibits 
behavior similar to the gas profile, in that the column den- 
sities at all radii drop over time. While the high-redshift 
version of the MFCS model matches the O VI observational 
data relatively well, once it has reached low redshift (the 
redshift of the actual observational dataset), it has fallen 
almost an order of magnitude. This behavior of ionic col- 
umn densities to decline by an order of magnitude at all 
radii is somewhat generic across all of our species, although 
some of the line/simulation combinations have a more sub- 
tle time-dependence belying the complex nature of the Saha 
equation. 



curves of the resulting galaxies as a metric of success. We 
revisit this analysis with the current models in Figure [7] 
Consistent with Paper I, simulations V (canonical run), LF 
and MF produce galaxies with very peaked rotation curves 
(Wcirc ~ 600 km s _1 ), clearly suffering from the angular mo- 
mentum problem. The simulation including cooling suppres- 
sion (MFCS) was successful at reducing this peaked value 
significantly to v c i TC ~ 400 km s _1 . The high feedback run 
resulted in a system with a very flat rotation curve and 
about 20% less total mass than the other halos. While our 
MFCS and HF simulations may appear as successful mod- 
els (i.e. realistic) based on this metric alone, the preceding 
analysis comparing ionic column densities makes it clearer 
that they do not predict significantly different absorber dis- 
tributions from the lower feedback runs (and fail to match 
observations). This fact alone reveals that QSO absorption 
line observations are an important constraint for feedback 
models to match in addition to traditional measures of suc- 
cess (e.g. rotation curves, star formation rates, Tully-Fisher 
Relation, Kennicutt-Schmidt Relation). 



4 DISCUSSION 

For all of the simulations, we see the same general behavior 
over each of the ionic radial profiles. The low-energy species 
of H I, Mg II, Si II, Si III, Si IV, and C IV display a core 
column density level out to a characteristic radius, Tinner, 
beyond which they fall steeply to a floor of trace levels in- 
dicating that these ions tracing cool gas are mostly asso- 
ciated with the disk. As we increase the ionization energy 
of the species, the core column density diminishes and the 
floor column density level increases indicating warmer, more 
distributed gas being traced. This trend continues to the 
highest energy ions of O VI and O VII, which display hot, 
extended gas in mostly flat radial profiles. 

It is clear from our results that these simulations do 
not produce column densities distributions of various ionic 
species matching QSO absorption observations. In general, 
our models predict column densities which are consistently 
at lower values and less extended than the observational 
studies, save for the predictions of the hottest gas in O VII 
absorption. Furthermore, there is only a hint of a trend in 
our data suggesting a path to resolve this discrepancy (e.g. 
by increasing feedback efficiency). So why do our models fail 
to agree with observations? 

One obvious problem with our method is that we only 
consider column densities of material within 1 Mpc of our 
galaxy. In reality, quasar sightlines pass through a much 
larger volume, and it is possible that additional intervening 
absorbers with velocities similar to our galaxy exist, which 
would increase the column densities in subsequent radial 
projections. However, this effect is likely to be small, so we 
seek additional explanations for these column density dis- 
crepancies. 



3.4 Rotation Curves 

In Paper I, we analyzed the efficacy of different subgrid 
models (including feedback prescriptions) using the rotation 



4.1 Matching O VI 

In order to better understand the shortcomings of our mod- 
els, we attempt to reconstruct the conditions necessary to 
reproduce the observational datasets. In particular, we try 
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Figure 6. Radial profiles for gas, metallicity and O VI as they change over time. The left and right sides display the evolution of radial 
profile for the V and MFCS runs respectively. 



to match the O VI observations, since of all the observa- 
tional datasets compared, they have the smallest errors and 
are the least likely to be susceptible to inaccuracies from 
photoionization (see section 4.2 below). As noted previously, 



our radial profiles of O VI column density are 1-2 orders of 



magnitude lower than the observations in the r ~ 10 — 100 
kpc region. 

Broadly speaking, one could imagine two explanations 
for this deficit. First, there could simply be an insufficient 
amount of metal-enriched gas at the projected radii in ques- 
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Figure 8. Radial profiles of O VI column density for a single data output at z = 0.2. The left and right sides show the V and MFCS 
runs respectively. The O VI column density of the original model is displayed in black, whereas the maximum O VI column density when 
we leave temperature as a free parameter is displayed in blue. 



tion, or second, the gas could be there in sufficient quantities, 
but not in the correct ionization state, perhaps because the 
temperature is too high or too low. O VI is particularly sus- 
ceptible to the thermal state of the gas, as it only exists in 
sufficient quantities in a relatively narrow temperature range 
(assuming collisional ionization). To explore which of these 
two explanations is most consistent with our simulations, 



we take the radial profiles of physical quantities generated 
in Figure [5] and systematically vary either the density, tem- 
perature, or the metallicity at each projected radius in order 
to maximize the resulting O VI column densities, holding all 
other quantities constant. Under the assumption that colli- 
sional ionization of O VI dominates, we fix the radiation 



intensity (see Section 4.2 1 
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We test this out on a single data output at z = 0.2 
in the MFCS run although the result holds at other times. 
First, we find that for a fixed metallicity and temperature, 
changing the density has little effect - this might be some- 
what counterintuitive, but is expected because we are keep- 
ing the column density constant, and O VI collisional ion- 
ization depends very little on density. The metallicity has 
a much larger effect, but we find that there is no solar or 
sub-solar metallicity values which will bring our column den- 
sities up to the observed levels (we limit the metallicity to 
solar on plausibility grounds). On the other hand, modifying 
the temperature (i.e. by lowering it to 3 x 10 5 K for most 
radii) while holding everything else fixed succeeds in pro- 
ducing sufficient column densities to match (or even exceed) 
the observational data, as seen in the right panel of Figure 

m 

This suggests that our MFCS model possesses sufficient 
quantities of enriched gas at large radii, but that insufficient 
amounts of cooler gas exist (necessary for O VI). Note that 
this procedure is not physically realistic - we are not free to 
reduce the gas temperature everywhere (this would decrease 
the pressure and hence change the gas distribution); how- 
ever, one solution would be to introduce multi-phase gas, 
preserving the total pressure but increasing the amount of 
gas in the 3 x 10 5 K temperature range. This result appears 
to apply to the other feedback simulations as well, that is, 
the median temperature is generally too high in our CGM, 
leading to depleted concentrations of these ions in preference 
to higher ionization states, but modifying the temperature 
allows us to match the observations. 

Interestingly, the results for the Vanilla run (shown in 
the left panel of Figure [8| are quite different - fine-tuning 
the temperature in this model cannot bring it in line with 
observations, since the metallicity falls off far too quickly at 
r ~ 50 kpc to allow for the modified temperature to compen- 
sate for it, further confirming that feedback and outflows are 
necessary in order to produce an observationally consistent 
simulation. 



4.2 Limitations 

To better understand the limitations of our models, we look 
at uncertainties related to photo-ionization and collisional 
ionization. To help focus this discussion, we separate our 
chosen species based on their primary ionization source. 
Generally the 'high' ions (N V, O VI, O VII) are thought to 
be associated with collisional ionization, whereas the 'low' 
ions (Mg II, Si II, Si III, Si IV, C IV) can be the result 
of both collisional and photoionization (Danforth & Shull 



2008 ) . Oppenheimer & Dave ( 2009 ) suggest that many of the 



weaker O VI absorbers may be due to photoionization, how- 
ever even they agree that absorbers with equivalent widths 
greater than 100 mA (most of our observational compari- 
son dataset) are predominantly due to collisions. With that 
division made, we can discuss the the ionization sources sep- 
arately. 



4-2.1 Photoionization 

Our low-energy species are partially ionized by absorption of 
energetic photons, so any problems in our modeling of pho- 
toionization will affect their column density profiles. Our 



study calculates photoionization for a metagalactic ionizing 
background in the optically thin regime, but it does not 
account for local photoionizing sources, such as energetic 
stars and galaxies within the volume itself. Inclusion of these 
sources will increase the ionization state of the various atoms 
present, thus potentially increasing (or decreasing) the col- 
umn density of the ions we study. Its dominant effect will be 
felt close to the disk of the galaxy where the bulk of the star 
formation is occurring and the local radiation field is high. In 



a similar study, Shen et al. (20121 included local photoioniz- 



ing effec ts with the use of a g alactocentric STARBURST99 
modeQ ( |Leitherer et al.|l999) for a SFR of 20 M yr -1 and 
found that local radiation was only important in the volume 
r < 45 kpc. Given that our SFRs are typically in the range 
of 10 Mq yr -1 , we predict an even smaller region where local 
photoionizing sources are important. However, the discrep- 
ancy between observations and the models occurs primarily 
at large radius, so while the lack of local photoionization 
may have some effect on our profiles it seems unlikely to 
explain much of the observed differences in column density 
profiles for O VI, H I, Mg II and possibly C IV at distance 
well beyond ri nncr . 

Cool gas at high densities will act to shield itself from 
photoionizing radiation. In this paper, we have treated all 
gas in the optically thin regime for absorption of the ioniz- 
ing background; however, gas with column densities iVm > 
10 17 ' 2 cm" 2 should be partially or completely shielded from 



radiation more energetic than 13.6 eV (e.g., Faucher-Giguere 
|fc Keres|2011[ ) . This will have the effect of lowering the ion- 
ization state of the various species present, increasing the 
column density of the lowest ionization species preferentially. 
Examination of the H I radial profile reveals that column 
densities of this magnitude only occur in the near-disk re- 
gion of the galaxy (r < ri nncr ). Interestingly, our lack of 
self-shielding may in some cases be balanced by our lack 
of local photoionizing sources, since they both dominate in 
the disk and can have counteracting effects depending on 
the ion. Either way, neither seem likely to have significant 
repercussions at large radii, in the region of the greatest 
observational discrepancies. 

Like many simulation groups, we use a Haardt-Madau 
photoionizing background (Haardt & Madau 1996 updated 
to the 2005 version), but there have been some advances in 
constraining the nature of this background and its sources in 
the universe since this time. A few improved approximations 



of the background have appeared recently (e.g. Faucher- 
Giguere et al. 1120091 iHaardt & Madau||2012l). Employing a 



different background has the potential to modify the ion- 
ization balance, so as to change the column density of gas 
in these low-energy ions, although it is unlikely that it will 
elevate our results by the order of magnitude difference re- 
quired to align them with observations. 

Current work is ongoing to include the aforementioned 
physics into our simulations. However, while these effects 
might make some differences in the simulated profiles, the 
changes will probably be confined to the near-disk region, 
failing to bring our models substantially closer to agreement 
with observations. 



4 http:/ /www. stsci.edu/science/starburst99/docs/default. htm 
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4-2.2 Collisional Ionization Equilibrium 

The use of Cloudy in calculating the ionization state of 
our various species assumes collisional ionization equilibrium 
(CIE), which is an ideal condition requiring r coo i >> tcie 
( |Dopita fc Sutherland] |2003[ ). The CIE assumption holds 
in high-density, low-temperature regions where the cooling 
time is short, but it can break down in the low-density re- 
gions of the IGM (typically r > r vir ) ( |Cen fc Fan"g||2006| ). 
Cen & Fang ( 2006 ) investigated the differences in O VI and 
O VII absorber strengths with and without the assumption 
of CIE and found that CIE overpredicts equivalent widths 
for these two ions by a small factor (less than 50%). No- 
tably, properly accounting for non-equilibrium effects would 
decrease our column densities of O VI and O VII, pulling 
them farther away from the observational sample. 

The high-energy species associated with warm and hot 
gas are primarily ionized due to collisions with ambient elec- 
trons. The rate of these collisions increases with electron 
density and density of the medium, so that clumpy regions 
of high density should have a large amount of collisional 
ionization. Thus, resolution becomes a key element in or- 
der to assure one can resolve clumps sufficiently for proper 
modeling of their ionization state. 



4-2.3 Spatial Resolution 

To investigate the impact of resolution and its potential non- 
linear effects on the collisional ionization equilibrium in our 
species, we analyzed several outputs in an attempt to iden- 
tify the characteristic size of absorbers in various different 
ionic species. We found that clumps of material at moderate 
to high density in our fields came in a variety of sizes, rang- 
ing from over 10 kpc down to 425 comoving pc (our spatial 
resolution limit), but there was no preferred size scale of ab- 
sorbers. This continuum of sizes indicates that there may be 
more gas clumps forming on subgrid scales not resolved in 
this study. 

There is evidence that absorber sizes in reality are below 
the resolution used in this study. The size of the background 
sources, the quasar accretion disks, are typically ~ 0.01 pc in 
size ( Kochanek et al.|2007 l, allowing one to observationally 
probe very small-scale absorbing structures. Using multiply- 



imaged lensed quasars, Rauch et al. ( 1999 1 find that low ion 



absorbers can have a length scale < 20 pc, a result further 



confirmed by Simcoe et al. (20061. Photoionization mod- 



eling by Schaye et al. (20071 indicates that high ion ab- 



sorbers (e.g. C IV) are typically larger but often still quite 
small, ~ 100 pc. While not directly CGM related, |Faucher-] 
Giguere et al. (20121 demonstrate that radiative shocks can 



develop in AGN outflows generating compact clumps with 
sizes as small as ~ 0.01 pc. It is possible such a mechanism 
might produce similar clumps in galactic wind systems. Such 
clumps are potentially significant as a source of multiphase 
structure and ion absorption below the resolution of our sim- 
ulations. It is likely that increasing the spatial resolution 
of the simulation would better resolve some of these small- 
scale structures, increasing the overall column density of all 
low ion absorbers. This might additionally have an effect on 
high ion absorbers, as long as their ionization potential had 
a characteristic temperature less than the virial tempera- 



ture of the halo. This effect might reduce the discrepancies 
between these simulations and observations of the CGM. 

To further investigate this issue, we modified outputs 
from the simulations by resampling them at decreased max- 
imum resolution (0.8, 1.6, 3.2, 6.4 and 13. kpc maximum 
spatial resolution instead of the 0.4 kpc used in the simula- 
tions) prior to post-processing them for the atomic species. 
This step is done on all fields (mass- weighted for the temper- 
ature and metallicity fields) in order to determine whether or 
not our atomic post-processing step reaches resolution con- 
vergence. This downsampling will show us where the multi- 
phase substructure exists in the simulation, since smoothing 
two multiphase structures together (e.g. cold clumps in a hot 
medium) will result in gas with different ion populations. 

Figure [9] displays the outcome of this study for O VI 
and H I column densities of an MFCS output at z = 0.2. 
Neutral hydrogen, which is very cuspy in the high-resolution 
simulation (seen as the solid black line) is spread out over 
an increasingly large region as the spatial resolution dimin- 
ishes. The core is particularly sensitive to this effect, as the 
cold core is artificially mixed with the hotter surrounding 
medium, it raises the core mean temperature, leading to 
hydrogen ionization and less H I. On the other hand, the 
O VI column density actually increases in this region be- 
cause the artificial warm component favors O VI over other 
oxygen ionization states. In each case, we observe that our 
change in resolution in post-processing has very little effect 
on the CGM at large radius (except for a small increase for 
O VI) suggesting that most of our multiphase substructure 
exists within r ~ 50 kpc. By this measure, we are converg- 
ing on behavior outside of this radius; however, this result 
could change if the underlying simulation had higher reso- 
lution in the halo, or if it employed a different prescription 
for generating outflows. Alternative feedback prescriptions 
which do not utilize thermal energy deposition (e.g. radia- 
tion pressure from massive stars, kinetic winds from super- 
novae blastwaves, etc.) may be more successful at generat- 
ing multiphase outflows by not obliterating the cold, fragile 
medium. 



4.3 Comparison with Previous Work 

A number of recent studies have investigated the state of 
the CGM in hydrodynamic galaxy simulations both in the 
high-redshift and low-redshift universe. Many such studies 
focus on neutral hydrogen in the CGM of simulated galaxies 
at high z.|Altay et al.|(|2C)ll[) post-process an OWLS (Schaye 



et al. 2010[ ) simulation at z — 3 to include a self-shielding 
correction and the formation of molecular hydrogen, demon- 
strating that the inclusion of these features causes a flatten- 
ing and steepening of the H I column density distribution 
function at Nm ~ 10 18 cm" 2 and N m ~ 10 



spectively. Faucher- Giguere & Keres ( 2011 1 examine an SPH 



simulation of a galaxy at z = 2 — 4 by ray-tracing the ioniz- 
ing background to correct for self-shielding, discovering low 
covering factors of high-density H I in the CGM consistent 



with those we produce in our models. Barnes et al. (20111 



investigate how modifying the intensity of galactic winds af- 
fects the H I distribution and subsequent Lyman-a emission 
from such systems at z — 3. Their momentum-driven wind 
feedback prescription yields H I distributions similar to our 
own MFCS run at that redshift. 
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Figure 9. Radial profiles of O VI and H I column density as a function of downgraded spatial resolution for a single data output of 
MFCS at z = 0.2. 



Other researchers have explored using metal ions as pos- 
sible probes of the CGM by post-processing their simulations 
to include additional species, as we do in the current study. 



Kacprzak et al. ( 2010 1 added a Mg II field to an AMR galaxy 



simulation at z = 1 (similar to our MF run), finding a Mg 
Il-absorber size distribution and Mg II r inncr = 20 kpc con- 



sistent with our own work. Fumagalli et al. (2011 1 examined 
outputs at z ~ 2 — 3 from an AMR simulation like our MF 
model, processing them to include self-shielding due to dust, 
radiative transfer from local ionizing sources, and H I and 
Si II fields. In agreement with our results, they determined 
that the concentration of Si II absorbers is an order of mag- 
nitude lower than observed values, suggesting they too may 
lack significant multiphase structures in the CGM. A suite of 
simulations modifying the feedback prescription to include 
AGN and wind-driven feedback was performed by |Tescari| 
et al. (20111, demonstrating that C IV statistics were con- 



siderably more sensitive to feedback implementation than H 
I in the intergalactic medium (IGM) at z — 2—3. While their 
conclusion agrees with our results in the outer parts of our 
halos (r > 200 kpc), we find that both H I and C IV are sen- 
sitive to feedback prescription in the CGM (20 < r < 200 
kpc) as exhibited in Figure [5] Goerdt et al. (20121 post- 
process an AMR simulation of z ~ 2.5 galaxy output to 
include a self-shielding correction and H I, O I, C II, C IV, 
Si II, Si IV, Mg II, and Fe II, observing the distribution of 
these metals versus galactocentric radius using a feedback 
model similar to our MF simulation. In general our results 
at that redshift align well with their study, except our C 
IV column densities are an order of magnitude greater than 
theirs, which might be a due to their self-shielding correction 
suppressing higher levels of ionization in the CGM. Lastly, 
|Shen et al.| ( 2012[) analyze a z = 2.8 output from the Eris 
Simulation ( Guedes et al.||2011[ | akin to the MFCS model, 
including a spherically-symmetric correction for local pho- 
toionizing sources and adding H I, C II, Si II, Si IV, C IV, 



and O VI in post-processing. They show the radial depen- 
dence of O VI column density, generally agreeing (ours is 
flatter) with our MFCS results at z ~ 3 (see Figure [61 . 

At low redshifts, there have only been a few studies 
probing the nature of the CGM in numerical simulations. 



Van De Voort & Schaye (20111 demonstrated that inflow- 



ing, cold gas lacks significant metals in the CGM, consistent 
with our vanilla run presented in Figure [5] Turning to H I, 
Fernandez et al. (20121 looked at a z — galaxy (similar 



to our MF model) and examined the source of H I feeding 
its CGM. They showed H I column densities and covering 
factors matching our own study, exhibiting a steep falloff of 
iVm at r ~ 10 kpc but an extended floor of material beyond 
r v i r . Ford et al. (20121 analyzed a series of SPH simulation 



outputs at z = 0.25 to investigate the column density of var- 
ious metal lines in proximity to a large samples of galaxies 
of different characteristic mass. The results for their fiducial 
feedback model (utilizing explicit momentum-driven winds) 
on 10 12 Mq halos match well against our MFCS run (as 
seen in their Figures 2, 11, and 12), and the apparent dis- 
crepancies in Figures 9 and 10 appear to be due to the addi- 
tional normalization they employ for dealing with multiple 



halos. Stinson et al. (2011 1 produced a suite of SPH simula- 



tions run to z = including two feedback intensities. Their 
"low" feedback run was similar to our LF run but included 
cooling suppression, whereas their "high" feedback run (in- 
tensified in order to integrate the effects of radiation from 
massive, young stars) was ~ 3x stronger than our HF run 
while still including a form of cooling suppression. The re- 
sulting outputs were processed to include O VI data, and 
then radial profiles were produced for both H I and O VI. 
The column density profiles for the "low" feedback run are 
consistent with our LF simulation; however, the profiles of 
their "high" feedback run have column densities exceeding 
any of our simulations. These profiles align well with the ob- 
servational datasets in H I and O VI, suggesting that their 
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success may be a result of using extremely high feedback in- 
tensities (i.e. 10 51 ergs released per 7.5 Mq of stars formed). 
We see a possible trend in our models, potentially indicating 
that intensifying feedback could increase column densities of 
various species in the CGM region (notably O VI) but this 
prediction is not clear cut. 

There exists an additional physical effect that might 
account for some of the discrepancy between the |Stinson| 
et al. (20111 work and our own. Much of the gas in these 



simulated galaxies resides at the virial temperature of the 
halo TUr oc Mhaio/niaio- Our simulated halo is M200 = 1.1 X 
10 12 Mq whereas the halo used in their high feedback run 
has a mass M V - 1T — 7 x 10 11 Mq. This equates to roughly 
to a 35% decrease in the virial temperature of their halo 
relative to ours. As already noted, the O VI transition has a 
relatively small temperature window, centered around T ~ 
3 x 10 5 K. A modest shift in the ambient temperature of a 
halo could yield a significantly increase in the amount of O 
VI present, potentially explaining some of the disagreement 
in the results of these two studies. 



5 CONCLUSIONS 

We perform the first AMR study to probe the CGM at 
low redshifts by comparing simulations against a large host 
of observational probes and the first study to look at the 
redshift evolution of the simulated CGM. We carry out 
cosmological simulations of a Milky- Way-massed galactic 
halo using a variety of stellar feedback prescriptions, post- 
processing the data to compute the column density distribu- 
tions of a wide range of atomic species, spanning an ioniza- 
tion energy range from 13.6 eV to 739 eV to effectively probe 
different phases of the CGM. We median-average these dis- 
tributions over many orientations and times (from 2 = 0.1 
to z = 0.5) both in order to make a cleaner comparison with 
data, and to decrease the impact of satellites on our result- 
ing radial profiles. Lastly, we compare these results against 
a number of different low-z observational datasets, correlat- 
ing absorber strength and impact parameter to the nearest 
galaxy galaxies. Our main results are as follows. 

• It is difficult for any of our feedback models to repro- 
duce the column density distribution of some of the observed 
absorbers - in particular, we cannot match observations of 
Mg II beyond r > 50 kpc, and O VI at any radii. On the 
other hand, we find that the Mg II distribution at smaller 
radii (r < 50 kpc) can be reproduced in some of the models. 
In addition, the predicted O VII absorption (which probes 
hot gas) in simulations including moderate or strong feed- 
back is in agreement with current observations. 

• We find that some form of feedback is required to get 
significant amounts of metals into the CGM beyond ~ 50 
kpc, but the total amount of gas in this region of the CGM 
is not strongly affected by feedback. In models which include 
some minimal amount of feedback, we find sufficient metals 
in our simulated CGM to reproduce observations, but a lack 
of multi-phase gas means that not enough of the gas has 
sufficiently low temperatures to produce the observed ions. 
This is clearest in the case of O VI, where the observations 
generally do not match any of our models because of this 
mismatched phase (i.e. temperature). 



• Turning thermal feedback up to very high values (even 
without cooling suppression) is successful at producing a 
galaxy with a flat rotation curve; however, it still fails to 
match the observational datasets from QSO absorption line 
studies (notably O VI) demonstrating that QSO observa- 
tions provide additional constraints on subgrid models be- 
yond traditional metrics of success. 

• There is significant redshift evolution of the column 
density distribution for our atomic species. In general, the 
column densities for most species drop by almost an order of 
magnitude across all radii in evolving from 2 ~ 3 to 2 = 0. 



Generally, we find that QSO observations place tight 
constraints on feedback models. Perhaps not surprisingly, 
the CGM which acts as both stellar fuel and the medium 
into which supernovae explode, contains a great deal of in- 
formation regarding the specific nature of star formation 
and stellar feedback. Probing this medium through QSO 
absorption line studies can tap into this information and 
reveal clear ways of differentiating between subgrid mod- 
els. Significantly sampled observations in any of H I, Si II, 
Si III, Si IV, or C IV identifying ri nncr , the transition ra- 
dius where column densities plummet in the CGM for these 
low-ionization ions, should potentially favor a specific type 
and intensity of stellar feedback. Furthermore, even subgrid 
models which can meet other criteria for success (e.g. flat 
rotation curve, matching the Kennicutt-Schmidt Law) need 
not match QSO absorption line statistics, as is demonstrated 
in the case of our HF run producing a flat rotation curve but 
failing to reproduce the the Mg II, H I and O VI observa- 
tional datasets. Thus, comparisons against QSO absorption 
lines place a constraint on galaxy simulations orthogonal to 
standard comparison modes. 

Future work on this topic includes addressing some of 
the limitations considered in Section [4. 2| As was discussed, 
it appears that we lack significant multiphase substructure 
in the halos of our galaxies, which is likely responsible for 
our inability to produce sufficient column densities of certain 
atomic species required to match observations (e.g. H I, Mg 
II, O VI). The use of subgrid models which deposit stellar 
feedback explosions as thermal energy may be responsible 
for generating unrealistic hot outflows. Relying solely on ram 
pressure from hot winds to distribute feedback energy could 
result in the spurious destruction of any existing cold halo 
gas and suppress the creation of a multiphase medium in 
the CGM. Recent work on alternative feedback prescriptions 
may provide a solution. For instance, Murray et al.| ( 2011 \ 
emphasize that radiation pressure from energetic massive 
stars on the surrounding dust-entrained gas may provide a 
way to get cool gas out to radii where it is commonly ob- 
served. There have been a few promising demonstrations of 
radiative feedback (e.g. Kim et al.|2 011 2012 Hopkins et al 



2012 Wise et al. 2012 1. Likewise, subgrid models of feed- 



back where supernova energy is deposited as kinetic winds 
Springel & Hernquist||2003| lAgertz et al. 



(e.g. _ 

cosmic rays (e.g. Miniati 2001, Jubelgas et al. 



20121 or as 



2008 Vazza 



11 ■ ■ ■ ■ ■ 

2012) might also avoid destroying a fragile multiphase 

medium. Additionally, steps to ensure higher resolution in 

the galactic halo may result in smaller clumps forming there; 

however, one potentially must modify refinement criteria to 

achieve substantial resolution in the low-density halo. 
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